Tilted Physics: A Cosmologically Dipole-Modulated Sky 
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Physical constants and cosmological parameters could vary with position. On the largest scales 
such variations would manifest themselves as gradients across our Hubble volume, leading to dipole- 
modulation of the cosmic microwave anisotropics. This generically leads to a correlation between 
adjacent multipoles in the spherical harmonics expansion of the sky, a distinctive signal which should 
be searched for in future data sets. 
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Introduction. — Physical laws are described by a set 
of fundamental dimensionless constants. These in- 
clude the electromagnetic fine structure constant, a c = 
e 2 / (4ireohc) ~ 1/137, and the much weaker gravitational 
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fine structure constant, 



a G = Gml/(Hc) ~ 6 x 10 
(here m p is the proton mass; see e.g. [1] and references 
therein). The standard model of cosmology (SMC) (see 
e.g. [2]) has an additional set of adjustable parameters, 
some of which may derive from the first group. However, 
we cannot yet calculate these parameters from the known 
physical laws, and indeed it is not clear that we will ever 
be able to do so in practice. 

It is possible that any of these constants or parameters 
could vary in spacetime. Such a discovery would point to 
physics beyond the standard paradigms; for example, a 
spatial variation could indicate that separate patches of 
the Universe have different realizations of physics, pos- 
sibly related to the string landscape. It is important 
to constrain such variations using available cosmological 
data, since one might expect effects to show up on the 
largest observable scales. 

There has been a significant amount of activity inves- 
tigating the potential time variation of a c (see Ref. [3] 
for a review). Recently, studies of high redshift quasars 
have led to a claim of evidence for a spatial gradient 
(or "tilt") in a c [4]. Specifically, a dipole anisotropy of 
the form Aa e /a e = B(to — t(z)) cos 6 + m was fitted to 
quasar spectra, where B and m are constants, to — t(z) is 
the look-back time (as a function of redshift z) , and is 
the angle with respect to the dipole orientation. It was 
claimed that B > at ~4cr significance, although this 
result is yet to be tested independently 

An equally surprising claim made within the last few 
years has been that there is a very large-scale "cosmic 
flow" [5, 6]. Such a flow would be in conflict with predic- 
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tions of the SMC, but could be produced by a large-scale 
isocurvature mode. 

A gradient in a c or any other parameter could be 
observed (or constrained) at an early epoch using the 
CMB. For some of the parameters (e.g. a c ), the imprint 
on the CMB arises from a modification of recombina- 
tion physics, giving an epoch of last scattering which is 
anisotropic. For other parameters (e.g. the primordial 
amplitude), a gradient would be directly reflected in the 
properties of the primary CMB anisotropics. For another 
class of parameters (e.g. the density parameters), a gradi- 
ent would correspond to a long-wavclcngth isocurvature 
mode. In this Letter we describe how one would search 
for this effect using the higher order multipoles of the 
CMB covariance matrix. 

Several previous studies have discussed the breaking of 
statistical isotropy in the primordial fluctuations (e.g. [7]) 
and others have looked at the effects and observability of 
dipole or quadrupolc multiplicative factors applied to the 
temperature field on the sky (e.g. [8, 9]). Although these 
ideas are related to ours, here we are considering spa- 
tial variations in the underlying parameters themselves. 
In Ref. [10], the effect of spatial variation of cosmolog- 
ical parameters on the CMB was studied, although the 
covariance matrix was not considered. 

Dipole modulation of physics. — To make the problem 
tractable within standard perturbation theory, we as- 
sume that there is little spatial variation of each pa- 
rameter inside the Hubble radius at the time of re- 
combination, i rcc . This means that insofar as the lo- 
cal physics is concerned, the gradient can be consid- 
ered a first-order perturbation, and hence we can use 
the Friedmann-Lemaitre-Robertson- Walker background 
to evolve the matter perturbations as usual up to t rcc , 
with the parameter value(s) depending on the location 
on the last scattering surface (LSS). Today, however, we 
are observing photons free streaming from much greater 
distances, so the spatial gradient may become noticable 
over these scales (see Fig. 1). 

The spatial variation could be of arbitrary complexity, 
but here we restrict our attention to a linear gradient 
across our Hubble volume. This is motivated partly by 
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to an "intrinsic" dipolc on the sky. On the other hand, 
modulation of other parameters, such as A Sl will clearly 
not correspond to isocurvature modes. When an isocur- 
vature mode is present, the evolution of the mode should 
properly be examined to late times, to capture any line- 
of-sight effects that our approach would miss. Therefore, 
we present the higher order correlations for /b with the 
caveat that the connection with isocurvature models is 
yet to be fully explored. 

CMB covariance matrix. — The CMB temperature 
anisotropics observed today are described by T(n) = 
T [1 + $(n)], where the T is the monopole. The fluc- 
tuations are typically expanded in terms of the spherical 
harmonics Y^ m (n) as 



FIG. 1: Spatial gradient in a cosmological parameter (repre- 
sented by the colour gradient), which is locally negligible at 
last scattering, but important over the free-streaming length 
scale. The filled circle shows (to scale) the comoving horizon 
size at last scattering. 



the suggestions of a gradient in a c and by the flow mea- 
surements, but can also be justified in that any residual 
inhomogeneity from before inflation would be expected 
to have been stretched until it was essentially linear in- 
side our Hubble volume. Therefore we characterize the 
anisotropy in parameter X at t lec by 



X(n) =X Q (1 + AX/X ) cos( 



(1) 



where n is the direction to the LSS and X$ is the 
monopole value. We find that such a dipole induces cor- 
relations between the £ and I + 1 multipole modes of the 
CMB covariance matrix, with an amplitude and spec- 
trum depending on X. This is a generic feature of a 
dipolc anisotropy — our motion with respect to the CMB 
rest frame [11, 12] or a Hubblc-scalc homogeneous mag- 
netic field [13, 14] also induce similar correlations. 

In the following we fix the monopoles Xo to their re- 
spective WMAP 7- year best-fit values [15]. As well as 
a , we consider the baryon matter fraction /b = ^b/^m 
at the LSS, and the amplitude and spectral index of the 
primordial fluctuations, A s and n s . A wider list could 
be considered, such as: m p /m G ; other particle couplings, 
mass ratios, and mixing angles; helium abundance; or ef- 
fective number of neutrino species. However, our method 
assumes that the gradient affects only the generation of 
primary anisotropics (i.e. anisotropics generated at last 
scattering), so it should be generalized for parameters 
such as the redshift to reionization or tensor contribu- 
tion, which effect anisotropics along the line of sight. 

Isocurvature mode. — Some variations may lead to an 
additional £ = 1 mode. In the case of motion with respect 
to the CMB frame, any difference between the dipole de- 
duced from aberration and the directly measured CMB 
dipole can be ascribed to an isocurvature mode. In the 
same way, gradients in parameters like /b will correspond 
to gradients in the relative densities of matter compo- 
nents, i.e. to isocurvature modes, and hence will give rise 



i>(n) = ^2 ^2 a imYe m (n) . 



(2) 



In the approximation that the anisotropies are primary, 
we can Taylor expand the temperature field in terms of 
the unmodulated (statistically isotropic) field $'(n) and 
parameter X to give 

$(n) = $'(n) + AX-^i^/ -Y w {n) , (3) 

where we have oriented the dipole along the polar direc- 
tion. The harmonic modes are then given by 



dX 



(4) 



with the coupling coefficients 



1 (£+l) 2 -m 2 
(2£+l)(2£ + Z) 

1 {£ 2 - m 2 ) 



(5) 



/ (2£- 1)(2£ + 1) ' 
From this one can compute the correlation matrix 
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to order AX, where {a\* m a\, m ,) = C e S u ,S mml . The only 
non-isotropic terms are the off-diagonal £,£±1 modes — 
as expected, a dipolc modulation couples multipolcs with 
\A£\ = 1. The derivative power spectra can easily be ex- 
tracted from standard CMB codes; we used CAMB [16], 
using the unlensed spectra. 

The off-diagonal terms due to the aberration caused 
by our motion with respect to the CMB rest frame have 
a similar form, also coupling modes with |A^| = 1. These 
are given by (see e.g. [11]) 



Ce^m.i+l.m ~ P{£ + 1) [pl+l ~ C g ) fu,m,£ 



+ l,m ? 



(7) 
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FIG. 2: Derivative power spectra AX dC'i/dX with 
AX/X = 1CT 3 for X = a c (solid, black), f h (dotted, red), 
A s (dashed, blue), and n s (dot-dash, green). 
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FIG. 3: Fractional error on AX/Xq for a cosmic variance 
limited experiment out to I = f max . Labeling is the same as 
in Fig. 2. 



where (3 = v/c « 1.23 x 10 3 . This is essentially 
a derivative of the Cg with respect to £ itself, since 



C, 



C e m dCi/dt. 



Variable a c example. — To choose a concrete example, 
let us focus on a c . We recap the basic recombination 
equation, considering the 3-level hydrogen atom. A full 
explanation of the symbols can be found in Ref. [17]. The 
rate equation is 



din z 



x c x p n H « H - 0h(1 - x v )c- hv ~' kT ™\ (8) 

1 + -R:hAh^h(1 - x p ) 
X 1 + Kk(Xk + l3 H )n H (l - xp) ' 

where x p = n p /riH is the ionization fraction, /?h the ion- 
ization rate, Ah the 2-photon rate, Kn = X%/(8itH(z)) 
the rcdshift factor, and an the recombination rate. 

To adapt the RECFAST recombination code, we con- 
vert all energies and rates according to their a c de- 
pendence. The modifications are as follows: A Q — > 
A Q (l + e)- 2 ; K n K H (l + e)' 6 ; a H a H (l + e) 5 ; 

->■ /3h (l + e) 5 e~ e2 ; and A H -> A H (l + e) 8 , where 
e = Aa e /a e cos 0. We found the changes to the ionization 
history to be consistent with other work. An increase in 
a e leads to faster recombination, which is primarily due 
to the decrease in the ionization rate. 

Results. — In Fig. 2 we show the derivative power spec- 
tra for several cosmological parameters with amplitude 
AX/Xo = 1CP 3 . The Ct spectrum has a stronger relative 
dependence on a c than other parameters, giving larger 
off-diagonal correlations. We can see that the spectra 
have different shapes, and hence the dipole modulation 
of various parameters could in principle be distinguished. 

Note that the spectrum for X = n s depends on the 
chosen value of the pivot scale, ko. This is because, when 
A s is held fixed, dCi/dn s will be small near fc . Different 
ko would correspond to physically different spectra; our 
calculations are for the case of ko = 0.05 Mpc -1 , which is 



representative of pivot scales corresponding to the high-£ 
region. Also, defining a gradient for n s might be difficult 
for scales approaching the Hubble scale today. However, 
most of the constraining power comes from much smaller 
scales, where such ambiguity does not arise. 

We now turn to the question of the error on A X/ Xq for 
a single parameter from an idealized (i.e. cosmic variance 
limited) experiment. This could be studied in more detail 
using a Fisher matrix formalism, but since it is unclear 
which set of parameter variations should be considered, 
we leave this study for the future. 

The correlation matrix C^^ mj £^-x,m 

has 2£ + 1 off- 
diagonal modes for each £, which are diagonal in m for 
a dipole orientated along the polar direction. The vari- 
ance of each mode is, to lowest order in AX, CgCt+i. 
One can therefore construct the total signal-to-noise of 
the correlations by summing over £ and m (see e.g. [12]) 



E E 

£=2 m=- 



r 2 



(9) 



+i 



For a general dipole direction, the harmonic modes 
are rotated according to a' tm = J2 m > ^l'm(^'/'.7) 8 <m', 
where <f>, ip, 7 are the Euler angles and D l m , m the Wigner 
D-matrix. Since we orient the (unrotated) dipole in the 
polar direction, 7 = 0. The correlations arc then redis- 
tributed between non-diagonal m modes, depending on 
the Euler angles. However, one can show that the S/N is 
invariant under rotations, so (9) holds for any direction. 

An estimate of the error on AX/Xq can be derived 
from the total S/N by a(AX/X a ) = AX/(X a S/N) [12]. 
In Fig. 3 we plot this error as a function of £ max . At 
4iax = 2000, which is roughly the Planck satellite cos- 
mic variance limit, a(AX/X ) w 10~ 4 for a c and ~10~ 3 
for the other cosmological parameters. This compares to 
cr(/3) - 2 x 10~ 4 for the aberration effect [11, 12]. 

Returning to the case of a e , if we assume that the 
dipole variation claimed from quasar spectra at z ~ 
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2 continues as a linear slope in look-back time, then 
Aa /a c sw 1.5 x at i rec . Supposing instead that 

the gradient is linear in comoving distance r(z), then 
Aa e /a e would be a factor of roughly two higher at the 
LSS. However, this is still below the detectable limit from 
measurements of the CMB according to our results. 

In practice, the variation among different parameter 
gradients may contain partial degeneracies, and hence 
disentangling various possibilities may be challenging. 
Also, for small AX, it may be difficult to distinguish the 
correlations from those due to the aberration effect. On 
the other hand, a parameter gradient producing correla- 
tions significantly larger than those of aberration should 
be easy to detect. Also, weaker gradients may be de- 
tectable in the future using the extra information con- 
tained in 21 cm surveys [18]. 

One can also ask what happens to CMB polarization 
anisotropics in the presence of gradients. This has been 
studied for the case of a dipole modulation in the 3- 
dimensional potential (e.g. [19]), which corresponds to 
our case X = A s , and one finds similar correlations be- 
tween neighbouring multipolcs in polarization. This ad- 
ditional data can improve sensitivity by a factor of order 
two for the abbcration effect [12], and we expect the same 
to hold true for the cosmological parameter variations. 

The case X = A s has recently been investigated in the 
WMAP data [9], The authors found marginal evidence 
for a dipole modulation with amplitude AA S /A S = 0.07 
at i < 60, but found the effect decreased at higher £. It 



would therefore be interesting to check the WMAP data 
for other parameters with different correlation spectra 
than the primordial amplitude. 

Conclusions. — The spatial variation of fundamental 
constants has been considered before (e.g., [20-22]), but 
here we have shown how a gradient across our Hubble 
patch generically gives rise to correlations between neigh- 
boring multipolcs in the CMB anisotropy spectrum. 

Obviously we could extend the dipole to the case 
of quadrupolar modulation, giving correlations between 
multipole moments with |A£| = 2, and in principle go 
beyond the quadrupole. However, the physical interpre- 
tation then becomes less clear, and one would need an 
underlying model to relate all the degrees of freedom. 

Returning to the simplest dipole case, if such a "tilt" 
was found in a parameter, then an explanation would 
have to be sought in the very early Universe. Inflation- 
ary models have been suggested which lead to anisotropic 
initial power spectra [23, 24] and specific models might 
also explain spatial gradients in parameters. Other pos- 
sibilities include modulated reheating (e.g. [25]), bubble 
collisions (e.g. [26]) and topological defects (e.g. [27]). 
Certainly it is worth testing future data sets for these 
correlations, since they could be a smoking gun for vari- 
ations in physics on the largest accessible scales. 

Acknowledgments. — This research was supported by 
the Natural Sciences and Engineering Research Council 
of Canada. We thank Ted Bunn, Arthur Kosowksy, Ali 
Narimani, and Kris Sigurdson for useful discussions. 



[1] A. Moss, A. Narimani, and D. Scott (2010), 

arXiv:1004.2066 [astro-ph.CO]. 
[2] D. Scott, Can. J. Phys. 84, 419 (2006), 

arXiv:astro-ph/0510731. 
[3] J.-P. Uzan (2010), arXiv: 1009.5514 [astro-ph.CO]. 
[4] J. K. Webb et al. (2010), arXiv: 1008.3907 [astro-ph.CO]. 
[5] A. Kashlinsky, F. Atrio-Barandela, D. Kocevski, 

and H. Ebeling, Astrophys. J. 686, L49 (2008), 

arXiv:0809.3734. 
[6] R. Watkins, H. A. Feldman, and M. J. Hudson, MNRAS 

392, 743 (2009), arXiv:0809.4041. 
[7] C. Gordon, W. Hu, D. Huterer, and T. M. Crawford, 

Phys. Rev. D72, 103002 (2005), arXiv:astro-ph/0509301. 
[8] H. Zheng and E. F. Bunn, Phys. Rev. D82, 063533 

(2010), arXiv:1003.5548 [astro-ph.CO]. 
[9] D. Hanson and A. Lewis, Phys. Rev. D80, 063004 (2009), 

arXiv:0908.0963 [astro-ph.CO]. 
[10] E. P. Donoghue and J. F. Donoghue, Phys. Rev. D71, 

043002 (2005), arXiv:astro-ph/0411237. 
[11] A. Kosowsky and T. Kahniashvili (2010), 

arXiv: 1007.4539 [astro-ph.CO]. 
[12] L. Amendola et al. (2010), arXiv:1008.1183 [astro- 
ph.CO]. 

[13] R. Durrer, T. Kahniashvili, and A. Yates, Phys. Rev. 

D58, 123004 (1998), arXiv:astro-ph/9807089. 
[14] A. Bernui and W. S. Hipolito-Ricaldi, MNRAS 389, 1453 

(2008), arXiv:0807.1076. 



[15] E. Komatsu et al. (2010), arXiv: 1001.4538 [astro-ph.CO]. 
[16] A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 

538 , 473 (2000), arXiv:astro-ph/9911177. 
[17] S. Seager, D. D. Sasselov, and D. Scott, Astrophys. J. 

523, LI (1999), arXiv:astro-ph/9909275. 
[18] A. Lewis and A. Challinor, Phys. Rev. D76, 083005 

(2007), arXiv:astro-ph/0702600. 
[19] C. Dvorkin, H. V. Peiris, and W. Hu, Phys. Rev. D77, 

063008 (2008), arXiv:0711.2321 [astro-ph]. 
[20] J. D. Barrow and C. O'Toole, MNRAS 322, 585 (2001), 

arXiv:astro-ph/9904116. 
[21] J. F. Donoghue, Journal of High Energy Physics 3, 52 

(2003), arXiv:hep-ph/0101130. 
[22] K. Sigurdson, A. Kurylov, and M. Kamionkowski, Phys. 

Rev. D68, 103509 (2003), arXiv:astro-ph/0306372. 
[23] A. E. Gumrukcuoglu, C. R. Contaldi, and M. Peloso, 

JCAP 0711, 005 (2007), arXiv:0707.4179 [astro-ph]. 
[24] L. Ackerman, S. M. Carroll, and M. B. Wise, Phys. Rev. 

D75, 083502 (2007), arXiv:astro-ph/0701357. 
[25] J. R. Bond, A. V. Frolov, Z. Huang, and L. Kofman, 

Phys. Rev. Lett. 103, 071301 (2009), arXiv:0903.3407 

[astro-ph.CO]. 

[26] S. Chang, M. Kleban, and T. S. Levi, JCAP 0904, 025 

(2009), arXiv:0810.5128 [hep-th]. 
[27] M. Cruz, E. Martmez-Gonzalez, P. Vielva, J. M. Diego, 

M. Hobson, and N. Turok, MNRAS 390, 913 (2008), 

arXiv:0804.2904. 



